Hopping thermoelectric transport in finite systems: boundary eff"ects 
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It is sliown tliat for tlie hopping regime, the thermopowers in both finite two-terminal and three- 
terminal systems are governed by the edges of the samples. This is due to the fact that the en- 
ergy transfer between a transport electron and a conducting terminal is determined by the site 
most strongly coupled to that terminal. One-dimensional systems with both nearest-neighbor and 
variable-range transport as well as certain types of two-dimensional systems, are considered. For a 
given sample, the changes in the thermopowers due to modifying the bulk are quite limited, com- 
pared with those of the conductance. When the small thermopower changes exist, their average 
over a large ensemble of mesoscopic samples will vanish. We also obtain the distribution of the 
thermopower in such an ensemble and show that its width approaches a finite limit with increas- 
ing sample length. This contrasts with the distribution of conductances in such systems, whose 
width vanishes in the long sample limit. Finally, we find that the thermal conductances in the 
three-terminal case have a boundary-dominated contribution, due to non-percolating conduction 
paths. This contribution can become dominant when the usual conductance is small enough. All 
our theoretical statements are backed by numerical computations. 

PACS numbers: 72.20.Ec, 72.20.Pa, 84.60.Rb 
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I. INTRODUCTION 

Achieving large thermopowers is a challenge to our 
understanding of electronic transport. At the same 
time, it is a crucial ingrcdientii^ for many energy- 
conversion devices. In good, wide band, conductors 
the thermopower, S, is typically very small, due to 
the approximate electron-hole symmetry. Breaking this 
symmetry is therefore important for obtaining large 
values of S. This happens in various poor or nar- 
row band^ conductors/ semiconductors, near the metal- 
insulator transition^"— and in the hopping regime i^ii^ 

Recent experiments addressed thermoelectric trans- 
port in mesoscopic systemSi ^^'^^ Besides their general 
interest, they may be relevant for small-scale thermo- 
electric devices j^ii^ii^ Especially in the hopping regime, 
where the electronic states are localized and discrete, 
electron-hole symmetry is usually broken in a given sam- 
ple, even if it is restored by averaging over many real- 
izations. This should result in relatively large, sample- 
specific, thermopower. In addition, the parasitic phonon 
heat conductivity can be reduced due to interfaces and 
sample shape and geometry effectsJ^ii^ 

Most of the studies on the thermoelectric effects in the 
hopping regime were devoted to bulk systems i^i^iiiS Re- 
cently, we discussed the thermoelectric transport prop- 
erties in finite one-dimensional (ID) systemsi^ where 
the boundary effect was found to be very important for 
the thermopower. In this work we follow up, substanti- 
ate and generalize that study. The importance of edge 
effects on the thermopowers will be highlighted. We 
shall consider both ID and 2D finite systems, which 



can be arbitrarily large. Since the conduction electrons 
have to exchange energy with a reservoir, the "three- 
terminal geometry"—"— naturally appears. In addition 
to the two electronic terminals which exchange charge 
and energy /heat, the third terminal is purely thermal and 
mainly exchanges energy with the conduction electrons. 
The three-terminal setup for ID finite systems is shown 
in Fig. [TJ The system, bridging two electronic terminals 
(leads) consists of a number of localized states (LSs) with 
random energies. The system is connected to the leads 
by (dominantly elastic) tunnel couplings. Electronic con- 
duction through the system is achieved via tunneling and 
phonon-assisted hopping. The setup can be realized, e.g., 
when the two electronic leads are suspended and the sys- 
tem is mounted on a (boson bath) substrate The com- 
plete description of the thermoelectric transport in the 
linear-response regime is a 3 x 3 transport matrix relating 
the three currents to three "forces" (or "affinities")^"— 
see the next section. 

In this paper we consider a noninteracting localized 
system with constant density of states and localiza- 
tion length, in an energy window {—Ec,Ec)^ We start 
with a short review of the three-terminal thermoelec- 
tric transport (Sec. |ll|). We show that the boundary ef- 
fect dominates the thermopower in nearest-neighbor hop- 
ping (NNH) ID systems (Sec. IIIip . The simple under- 
lying physics is illustrated via the solution of a simple 
3-site model (Sec. IIII Al and Appendix [C]). Then, longer 
NNHID systems are considered (Sec. IIITB|) . In Sec. |lVl 
simplified types of 2D systems are treated. Because of 
incomplete averaging at the boundaries there can be a 
finite and fluctuating thermopower even for a very large 
system. The situation for variable-range hopping (VRH) 
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FIG. 1. (Color online) Schematic of the three-terminal geom- 
etry. A ID system (blue region) consisting of a number of LSs 
is connected to three terminals: the left and right ones (green) 
are electronic leads which have their own temperatures and 
chemical potentials, the terminal below (brown) is a thermal 
terminal (a phonon bath in this work) at temperature Tp. 



is discussed in Sec. |Vl In Sec. IV Al we shall first fo- 
cus on the ID case and then consider smaU 2D systems, 
whose width and length are comparable. In these cases 
the bulk effect plays a limited role, which will tend to 
disappear when the system size increases. When the sys- 
tem size is increased further, one would expect both the 
thermopower and its fluctuations to eventually go to zero 
in the macroscopic limit as a consequence of the particle- 
hole symmetry being restored with averaging. We find 
however that the fluctuations remain finite for an arbi- 
trarily long system, as long as its transverse size is finite. 
Finally, we discuss special effects of "non-spanning elec- 
tronic paths" (that do not transport charge) to the ther- 
mal conductances (Sec. IV Bp . Our statements are backed 
up numerically. The numerical scheme is explained in 
Appendix \^ and the contributions of various conduct- 
ing paths are compared in Appendix [BJ Finally, Sec. IVII 
includes a short summary and conclusions. 



II. THERMOELECTRIC TRANSPORT 
THROUGH LOCALIZED STATES IN THE 
THREE-TERMINAL GEOMETRY 

The study of hopping thermoelectric transport in 
three-terminal geometry was done in Refs. [l^ and ITtI . 
For completeness we summarize here the basic formula- 
tion of the problem. The hopping transitions between the 
LSs are assisted by phonons. We focus on the situation 
when the on-site Coulomb interaction is so strong that 
each LS can only be occupied by at most one electron. 
The inter-site Coulomb interaction may lead to interest- 
ing effects but these will not be discussed here»^ The 
hopping rate from LS i to LS j for, say, Sj > e^, is given 
by the Fermi golden-rule as 

r,^, = 271^ |a„q|2^(e, - - ^q)/,(l - /,)iY„-, 
q 

f,)m, . (1) 



Here, a^q is the electron-phonon interaction matrix el- 
ement between the two LSs, £j and Si are the energy 
of the LSs i and respectively, Wq is the phonon en- 
ergy, fj and fi are the occupation probability on the 
LSs j and i, respectively, and Nij is the phonon dis- 
tribution at the energy Wq ~ \ej ^ ei\. The phonon 
distribution in the system is determined by the phonon 
bath, Nij = [cxp(|ej - ei\/{kBTp)) - 1]"^ At large dis- 
tances the overlap of the wavefunctions of the two LSs 
is exponentially small. The asymptotic behavior of ^ij 
is thus ~ 7ep exp(— 2ry/^) where = jr^ — r^j is 
the distance between the LSs with Yj and denoting 
their position vectors, ^ is the localization length, and 
7ep is proportional to the electron-phonon coupling and 
the phonon density of states. The hopping rate from LS 
i to the left lead is 

y^^L^l^LI^{l- Ids,)) , (2) 

where = 27r|Q;iLpPL with standing for the cou- 
pling between the LS i and the extended states in the left 
lead of which pi is the density of states. We focus on the 
situation where 7^^ does not depend on the energy Si (i.e., 
no particle-hole asymmetry). Jl stands for the distribu- 
tion in the left lead, fhis) = [exp(e — yUi)/(fcsTL)) + l]^"^. 
The transition rate from i to the right lead can be writ- 
ten down similarly. The asymptotic behavior of 7^^ 
and 7iij at large distances is also exponential, 7iL ^ 
7eexp(-2riL/^) [7^^ ~ 7e exp(-2rjfl/^)] where r^i [rm] 
is the distance between LS i and the left (right) lead and 
7e scales with the tunnel coupling strength between LSs 
and the leads. The electric current flowing from i to j is 

li^j = e{T,^j - Tj^,) , (3) 

with e being the charge of the carrier. The electric cur- 
rent flowing from i to the left (right) lead is calculated 
similarly. At steady state, according to Kirchhoff's cur- 
rent law, 

h^j + h^L + h^R = 0. (4) 



which is also a statement of probability conservation. 
The steady-state distribution on each LS is obtained by 
solving Eqs. H]), (0), ©, and Q. 

Before formulating the currents among the three ter- 
minals in terms of the transition rates, we present a 
thermodynamic analysis in the linear-response regime. 
There are three heat currents flowing into each termi- 
nal, Ql, Qfi, and Qp, as well as two particle cur- 
rents flowing into the electronic terminals, and Nr. 
The heat currents are related to the energy and par- 
ticle currents according to the thermodynamic relation 
Qi ~ Ei — HiNi for i = L,R and Qp = Ep for the 
phononic terminal with Ei being the energy current flow- 
ing into terminal i. Particle and energy conservation ren- 
ders X]i=L R^i Ei=L R p ^ 0- Hence there are 
only three independent currents which are the electric 
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current If. = cNr = —cNl and two heat currents. For 
the latter one can choose Ql and Q/j, or any two hnear 
independent combinations of them. We shall adopt the 
convention introduced in our previous worki^ and choose 
the following heat currents^iii^ 

ih = \iQK -QL), i^Q ^-Qp = Ql + Qr. (5) 

In the linear-response regime the entropy production rate 
S is given by TS = T[{Ql/Tl) + {Qr/Th) + Qp/Tp)] = 
Ie{Sfi/e)+I^{S/T)+I^''{AT/T), where T is the common 
(equilibrium) temperature of the setup. This relation 
identifies the three "forces" (affinities) conjugated to the 
three currents, = fiL — Mi?, ST = — Tp, and AT = 
Tp — {TL+Tfj) /2. The phenomenological linear-transport 
equation which satisfies the Onsager reciprocity relations 
is theni^ 



5T/T 



L2 L3 KpJ \aT/T J 



(6) 



In the three-terminal geometry, besides the normal ther- 
mopower S = Li/{TG) there is the three-terminal ther- 
mopower Sp = L2/{TG) which converts the temperature 
difference AT to voltage (and vice versa) j ^^i^^ 

We now formulate the currents /g, /g, and in terms 
of microscopic quantities. The electric current is given by 



(7) 



The heat currents Iq and can be obtained from Ql 
and Qr, 



where ^ is the equilibrium chemical potential. 

In the linear-response regime the current between two 
LSs can be written as 



h^,^G,,{U,-Uj±U,, 



(9) 



where the signs -I- and — are for Sj > Si and Sj < Si, 
respectively, and the conductance of the bond (ij) is 
Gi, = (eV(fcBT))7,,/0(l - /j')(iVO +1/2 t 1/2) (the su- 
perscript is used to denote an equilibrium distribution 
function). InEq. ® C/. = {kpT/eM- f^)/{f^{l- f^)) 
and = {kBT/e){N,.j - N^^) / {N^j{N^^ + 1)). To com- 
plete the description of the resistor network we also write 
down the current between the leads and the LSs 



h^L = e(r,^L - Ti^O = G,L [U., - UL{e^)] 



(10) 



where G,l = {e" / {kBT))^^Lff{l ~ /°) and C/i(e,) = 
{kBT/e)[fL{e,) - /0]/(/0(l - /°)).^ We shall adopt the 
widely-used approximation, valid in the rather broad 



regimes discussed in Refs. |22h24 which yields 



2rij ~ + \ej - mI + \ei ~ Sj 



Gij ~ Go exp 1^ - 

GiL ^ Go exp ( - 



e 2kBT 



(11) 



where Go ^ e^lep/ {ksT) ^ e^^eil [kBT) has been in- 
troduced. This coefficient sets the scale of the whole 
conductance and will not play a role in the subsequent 
discussions. The resistor network described above has 
been investigated for the case AT = a long time agoi^ 
We have recently considered the effect of the term C/y , 
which arises due to a finite small AT>i^. A scheme for the 
numerical solution of the resistor network is presented in 
Appendix The above formalism reduces to the Miller- 
Abrahams resistor network mode l^^i^^ when there is no 
temperature difference. 



III. BOUNDARY EFFECT IN ID NNH 
SYSTEMS 

A. A simple three-site ID NNH system 

To demonstrate the boundary effect in the NNH regime 
wc study a simple model system which consists of just 
three LSs. Consider the situation where LSs 1 and 3 
are strongly coupled, by elastic transitions, to the lead 
continua on the left and on the right, respectively so 
that /i = /l(£i) and /a = /i?,(e3). That is, the con- 
ductances Gil and Gsij are much larger than the other 
conductances. It is also assumed that the tunneling con- 
ductances between LS2 and the leads are so small that 
the transport through the system is dominated by the 
hopping path illustrated in Fig. [2][a). The condition for 
this to be true is analyzed in detail in Appendix |B] This 
model system can be realized in experiments by, say, a 
serially coupled three-quantum-dots 1^ For concreteness, 
consider the situation when ei,e3 < £2- U2 is then de- 
termined by Ii^2 = ^2^-3 where 



/l^2 =Gi2(C/l -C/2 + C/12) , 
-^2^3 = G23(?72 — C^a — U23) , 

and consequently 

Gl2(C/l+C/l2) + G23(C/3+[/23) 



and 



U2 



G12 + G; 



23 



(12) 



(13) 



G12G23 
G12 + G23 



■[C/3 + t/23-C/l-C/i2] , (14) 



where le — Ii^2 is the total electric current. Express- 
ing the J7's as functions of the chemical potential and 
temperature differences gives 

J G12G23 r^/i , £31 <5T W31 AT 

-^e „ „ [— + — — + —\, 15) 

G12 + G23 e el el 
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(a) 



(b) 
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FIG. 2. (Color online) ID NNH. (a): Schematic of the three- 
site model in the three-terminal geometry. The three red dots 
stand for the three LSs. The arrows represent the hopping 
transitions in the dominant hopping path. The vertical (hor- 
izontal) direction stands for the energy (position). The color 
densities in the left and right sides represent the electronic 
distribution in the left and right leads, respectively. In the 
figure the left lead has higher chemical potential and temper- 
ature, (b) The effect of decreasing the conductances in the 
middle of the system by a factor Tg on the conductance and 
thermopower of a ID NNH system. X = G or S represents 
the conductance or the thermopower, X' / X denotes the ratio 
of the conductance or thermopower after modifying the mid- 
dle part of the system to its original value. Solid (dashed) 
curves are for X = S [X — G) for a three-site model. The • 
(A) points denote X = 5* (X = G) for a ID NNH model with 
31 LSs. Note that for both models S' /S = 1. The parame- 
ters for the 31 LSs ID NNH model are; ksT =20, n = 0, and 
^ = 0.001. The LSs are located at the sites of a ID lattice with 
a periodicity of 0.008. The energy is uniformly distributed in 
the energy window {—Ec,Ec) with Ec = 60. The Mott hop- 
ping length is .Jifpk^ 12 = 0.004. This value, half of the 
nearest-neighbor distance, makes all the hopping processes 
except the NN one sufficiently small. By "decreasing the con- 
ductances in the middle of the system" we mean decreasing 
the conductances of the connections among the 1st, 2nd, and 
3rd LSs, while for the 31 LSs ID NNH model it means de- 
creasing the conductance of all bonds between the 11th and 
21st LSs by a factor of Vg. 



where we have denoted £31 = (ei + £3)/2 — /i and wai 
£3 — £1. Using Eqs. ([7]) and ([5]) one finds 



Inserting these results into Eq. ([6]) yields the transport 
coefficients in the linear-response regime, 



ii = G ( ^ ) , L2 



G 



with 



G 



G12G2 



Gi 



G2 



(17) 



(18) 



and confirms the Onsager reciprocity relations. In Ap- 
pendix [C] we reproduce these results by a probabilistic 
analysis. 

Remarkably, the thermopower S = Li/iTG) as well as 
the coefficients L2/G, L^/G, K^/G, and ifpe/G depend 
all only on £1 and £3, i.e., the energies of the LSs at the 
boundaries. The site energy of the central level does not 
affect these quantities. The thermoelectric properties are 
completely determined at the boundaries. On the con- 
trary, the bulk (in this simple example, the central level) 
can affect the conductance of the system. Clearly, we can 
replace the middle site by a more complicated construc- 
tion. As long as it is coupled to the boundary sites in 
the same way, the properties of this mid-system do not 
matter for the above-mentioned transport coefficients! 

To illustrate the 3-site case, we have numerically com- 
puted the conductance G and the thermopower S* of a 
three-site system as a function of the decrease of the con- 
ductance of the bonds (1,2) and (2,3). Namely, we have 
determined the conductance and the thermopower of the 
system when G12 — ?• Gi2/rg and G23 — > G2zlrg. We plot 
the ratio of the new conductance G' to the original one G 
as well as the ratio S' / S as functions of the scale factor rg 
in Fig.[2ljb). It is seen that although the conductance de- 
creases significantly with increasing rg, the thermopower 
remains unchanged, S' / S = 1. More complicated models 
will be discussed below. 



B. Longer ID NNH systems 

We now extend the discussion to longer ID hopping 
systems. Nearest-neighbor hopping in a chain of LSs is 
accomplished via electron transits into the left (right) 
lead only through the leftmost (rightmost) LS, having 
energies £/ (£,■). Therefore from Eqs. ([7]) and ^ 



Ql 



Je, Vii — ; 



and one readily finds 



(19) 



(20) 



Ql 



£1 - M 



Qr 



£3 - M 



(16) 



with £,.f ~ {si + £r)/2 — ^1 and ujri = £r — £i- Inter- 
estingly enough, the thermoelectric properties can be de- 
duced without solving the resistor network. For example. 
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when 7^ and ST = AT = 0, one has le = G5fi/e 
and then by Eq. (^0)1 Iq ~ {ere/e)Ie = GeriS^JL/e^ and 
/q^ = {u}re./e)Ie = GuriS^I . Therefore 



Li = G 



is 



(21) 



Analyzing the situations when 5T ^ and = AT = 
and when AT ^ and 5^ = ST = and exploiting the 
Onsager reciprocity relations one obtains 



e 



pe 



(22) 



Again the thermoelectric properties are completely de- 
termined at the boundaries. This is clearly manifested in 
Fig. [2{b) for a NNH system with 31 LSs. It is seen that 
the thermopower is immune to the change in the middle 
part of the system while the conductance is considerably 
affected. Remarkably, this also implies that the ther- 
mopower (and other thermoelectric coefficients) is finite 
and random (as long as ei and Sr are finite and random) 
regardless of the bulk properties. This even persists, in 
general, for very long systems where the particle-hole 
asymmetry is negligible, as long as the edge sites have 
definite energies (and their sum does not happen to van- 
ish exactly). Therefore the particle- hole symmetry may 
no longer dictate a zero sample-specific thermopower in 
the macroscopic limit in ID NNH systems. 



IV. 



SIMPLIFIED 2D NNH SYSTEMS 



We start by studying a situation where the boundary 
effect fully dominates thermoelectric properties: When 
the electronic leads are geometrically sharp (as with a 
high resolution STM configuration) and each of them is 
coupled strongly only with a single LS as illustrated in 
Fig. |3l(a). Specifically the left lead is coupled with a LS 
of energy ee while the right one with a LS of energy Er- 
in this way the relation between the two heat currents 
(Iq and Iq^) and the electric current Ig is given again by 
Eqs. (|20p . Following the same logic as that applied for 
ID NNH systems one again obtains Eqs. ([2^ . Therefore 
the thermoelectric properties are completely determined 
by the boundary LSs (i.e., the LSs coupled strongly with 
the two electronic leads) in this case as well. 

Next we consider another special type of 2D systems 
that consist of parallel ID hopping chains where there 
is no transport between different chains as sketched in 
Fig. Efb). The heat and electric currents arc given by 



k k 



{k) 



(a) 



TrUr 









Tl 




Tr 




• • • • M 


I'r 






M • • •• • • • 





(b) 



phonon bath 



FIG. 3. (Color online) Three-terminal hopping transport in a 
2D NNH system, (a): Illustration of the situation when the 
two electronic terminals are geometrically sharp. Each ter- 
minal is strongly coupled only with a single LS with energy 
ee and Sr on the left and right, respectively. A possible hop- 
ping path is illustrated in the figure by the arrows, (b): A 
2D system consisting of a series of parallel hopping chains be- 
tween the two electronic leads. The hopping between different 
chains is negligible. 



where the superscript k denotes the fc-th chain. One then 
has 



Li=G 



G 



e 



- \ 2\ 



k2 




, (23) 

where G = "^j. Gk with Gk being the conductance of the 
fc-th chain and 



EkGk 



(24) 



Altering the central region will modify the conductances 
Gfc. Due to the random nature of the system this mod- 
ification varies among different chains and changes the 
averaging in Eq. (P^ . Although Sri and LUre for each 
chain do not change, the averaged values in Eqs. ((23)) 
do. Therefore the thermoelectric properties also depend 
on the bulk in this type of 2D systems. This is also true 
for other types of 2D NNH systems where the backbone 
consists of parallel hopping paths. Even when there are 
connections between those parallel paths, if the current 
distribution at the boundaries can be considerably af- 
fected by the bulk, the bulk effect cannot be ignored. 
However for sufficiently long 2D systems the current dis- 
tribution at the boundaries should not be affected by the 
far away bulk. In general there is a correlation length Leo 
in hopping systems beyond which the spatial current dis- 
tributions are uncorrelatedi^ In ID NNH systems L^o is 
on the order of the distance between the adjacent LSs. In 
2D NNH Leo can be much larger but still finite. In Fig.|3] 
we show how the conductance and the thermopower are 
affected by changes in the middle part of a long and 
a short 2D NNH system. In the long system the ther- 
mopower is almost unchanged whereas in the short one 
it is significantly modified (but, possibly much less than 
the conductance, which can be, for example, affected ex- 
ponentially). Leo for the chosen parameters is estimated 
as ^ 6 times the distance between the adjacent LSs. 
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(a) 



FIG. 4. (Color online) The effect of decreasing the conduc- 
tances in the middle of the system by a factor of rg on the 
conductance and thermopower, for a 2D NNH single realiza- 
tion. X = G or S represents the conductance or the ther- 
mopower, and X' /X denotes the ratio of the conductance 
or thermopower after changing the middle to the original 
value. The curves with A are for the thermopower while 
the ones with • are for the conductance. The parameters are: 
ksT =20, = 0, C = 0.01, = 20, and p = 40. The 
Mott hopping length is [S,/ {pkBT)Y'''^ /2 = 0.012. The LSs 
are put on a 2D lattice with periodicity 0.025. The num- 
ber of lattice sites perpendicular to the transport direction is 
Ny = 20. For the longer system (a) the number of sites along 
the transport direction is N^. = 31. For the shorter system (b) 
Nx ~ 15. Each lattice site is denoted by two indices (ix,iy) 
withij, e [-{Nx-l)/2,{Nx-l)/2] and iy £ [^,Ny]. The term 
"conductance in the middle" stands for the conductance be- 
tween two LSs {ixi,iyi) and {ix2,iy2) where ixi G [^5,5] and 
ix2 e [-5, 5]. 



We now discuss the macroscopic limit for these 2D 
hopping systems. For sufficiently long systems, the fluc- 
tuation in the conductance of each chain becomes very 
small.— The averages in Eqs. ((23|) become the averages 
over the energies at the boundaries. If the energies of 
LSs at the boundaries obey the same (sufficiently broad) 
distribution as the bulk ones, then 



Jdspje)... 
I dep{e) 



(25) 



with p{£) being the density of states of the involved LSs. 
When the density of states is symmetric with respect to 
the chemical potential, one finds 



Li = 0, L2 = 0, i3 = 0, 



1 (£2) 
-G— , Ape 



2G 



(26) 



with the average given by Eq. ([25)) . The particle- hole 
symmetry indeed leads to zero thermopowers in the bulk 
limit. In deriving Eqs. (f26| we have used {ei) = {Sr) 
and (e^) = (e^) = (e^), and have taken into account 
the fact that for very long ID systems there should be 
no correlation between ei and Er- When (e^) = (e^) then 
^2 = ^3 = in the bulk limit even when particle-hole 
symmetry is broken, such that Li ^ Q. In fact Sp = 
L2/{TG) and L3 have nothing to do with the particle- 
hole asymmetry but are related to the (Icft-to-right) in- 
version asymmetry in the sense that Sp oc (e^) — {si) and 
L,^{el)-{sl). 



V. BOUNDARY EFFECT IN ID VRH SYSTEMS 

Here we will not assume that only a single LS is 
strongly coupled to each lead. Quite generally, all LSs 
located within a distance from the lead smaller than or 
comparable to the Mott hopping distance Rm can be 
considerably coupled to that lead, with the conductance 
of the connection given by Eq. ([TT|). Other LSs, which 
are not coupled directly with the leads have much lower 
conductances, due to the exponential decay of Gil and 
Gm with the distance til and rm. This implies that the 
boundary effect is somewhat weakened. To study this 
situation, we consider sufficiently long ID VRH systems 
whose length L is much larger than Rm and denote by 
"boundaries" the regions that are within a distance of a 
few i?M's from the leads. We shall find that the bound- 
ary effect on the VRH thermopower is still important. 

Specifically for VRH systems the current flowing into 
each lead comes mainly from the LSs in the boundary 
regions. The summations in Eqs. ([7]) and ([5|) are then 
reduced to summations over those LSs. Accordingly, the 
thermopowers can be written as 



where 



{Si) 



Sp = ^ {{sr) - (ee)) 



{ee = 



(27) 



(28) 



with Y^'^ (E'Z) being restricted to the LSs in the right 
(left) boundary region. In ID NNH systems the summa- 
tion is restricted to a single LS that is coupled strongly 
with each lead, while in ID VRH systems there are more 
than one such LSs. Nevertheless, whenever the number 
of LSs involved in each summation is not too large the 
thermopowers will be finite and will fluctuate regardless 
of the particle-hole symmetry in the bulk. The boundary 
part of the backbone picture is drawn in Fig. [SJa). Simi- 
lar to NNH 2D systems, the weights of the various i's In 
Eqs. [28l(e,g, li-^L) does depend on the bulk. Therefore 
in a given sample, there will again exist some limited de- 
pendence of the thermopowers on the bulk. However, this 
will be averaged out in an ensemble of many realizations 
of the sample. 



A. Thermopowers 

We now turn to the thermoelectric properties of VRH 
systems in ID. The conductance of a ID system is mainly 
suppressed by the "breaks" rendering the Mott VRH 
picture not entirely applicable. At low temperatures the 
characteristic conductance of a ID VRH system of length 
2L is 



1.29 



GiD = GqC 



(29) 
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where 



1/2 



, TM = 2roln 



(30) 



with ksTo = (p^) ^ and p denoting the density of (local- 
ized) states, is the solution o&^iM 



To \ ^ J 



(31) 



The current mainly flows in the backbone of the resistor 
network which mostly consists of connections with con- 
ductance higher than or comparable with Gid^St— The 
typical hopping length and energy are Rm = ']C/2 and 
Em = riksT . Below we use the energy scale 



(32) 



which does not depend on the system length 2L. For 
example, T in units of this scale will be seen to be relevant 
for the thermopower fluctuations, see Fig. [7Kb). 

The boundary effect is detected by comparing the ther- 
mopower S (and other thermoelectric coefficients) of a 
random system and that of the same configuration but 
with the central region modified, S' . Specifically, we ap- 
ply the following modification: — > 10~^Gjj , if both 
i and j are in the central part. The boundary effect is 
monitored by | ln(|S"/S'|)|. Concomitantly we compute 
the conductances of the original and the modified sys- 
tems, G and G', and monitor the change via | ln(G'/G)|. 
If I ln(|S"/S'|)| is very small (i.e., S' is almost the same as 
S) then the thermopower is insensitive to the bulk and 
the boundary effect dominates. We model the localized 
electron system by a number of LSs located at random 
positions and and having energies which are uniformly 
distributed in the ranges {—L,L) and {—Ec,Ec), respec- 
tively. The central region is taken as x G {—Rm,Rm)- 
The linear-response transport coefficients are computed 
using the method described in Appendix [Al Note that for 
the numerics we use dimensionless energy and tempera- 
ture, with fcs = 1. For a given system, the appropriate 
energy unit can be introduced. 

The averages of the two quantities over 10^ random 
configurations are plotted in Fig. ^h) and (c). It 
is seen that the conductance is considerably modified, 
(I ln(G'/G)|) > 1. In contrast the change in thermopower 
is much smaller, especially when the distance between the 
central region and the boundary L—Rm exceeds the hop- 
ping length Rm- (| IndS'VS'DI) decays rapidly with the 
distance L — Rm and soon becomes negligible. 

For the choice of the central region adopted in the 
figure (from —Rm to Rm), the change of the conduc- 
tance is not dramatic, e.g., G' /G = 1/4 for the last point 
(N=120). But if the central region is taken to be from 
-2i?M to 2Rm then G'/G ^ 1/22. The relative changes 
in thermopower in the former and latter cases are, how- 
ever, no larger than 8 x 10~^ and 2 x 10~^ respectively. 

The three-terminal thermopower S'p, shown in Fig.[SKc, 
behaves similarly. Therefore the probability that the LSs 



far away from the boundaries can affect the average ther- 
mopowers is very small. This also indicates that the cor- 
relation length, Leo, giving the scale over which a local 
change in the network influences the conducting path, in 
ID VRH system is only a few hopping lengths. 

The question naturally arises, what happens in a spe- 
cific sample? We find that there again the boundary 
effect can be dominant. In Fig. ^ we plot the change of 
the thermopower and the conductance for two systems as 
a function of the increase of the resistance in the middle 
part of the sample. For the longer system it is seen that 
the thermopower is unaffected while the conductance de- 
creases by almost three orders of magnitude. This is the 
situation when the distance between the boundary and 
the middle exceeds the correlation length Leo so that the 
bulk affects the thermopower negligibly. However for a 
shorter system a change in the central part can affect 
both the thermopower and the conductance. Neverthe- 
less the change in the conductance is still much more 
significant than that in the thermopower. 

Interestingly enough, our analysis points out that the 
thermopower has unexpectedly large sample-to-sample 
fiuctuations even for very large samples. This is very dif- 
ferent from the vanishing of the VRH conductance fluctu- 
ations for increasing-length samples. No matter whether 
the bulk effect is important or not, as long as the number 
of the LSs involved in the summation in Eqs. (j28p is flnite, 
the thcrmopowers have a finite and fluctuating value due 
to insufficient averaging. For ID VRH, the LSs involved 
in the averaging, i.e., those with r^^ < Rm (oTrm < Rm) 
and \Ei\ < Em are typically just a fewi^ Thus the fluc- 
tuations of the thermopower can be rather large. To 
check this, we computed the variance of the thermopower 
as a function of the length of the system. The results 
are shown in Fig. [TJa). Indeed the variance of the ther- 
mopower remains considerably large and attains a con- 
stant value for very long systems. The variance of the 
thermopower in very long systems increases with increas- 
ing temperature (decreasing Rjm). The appearance of a 
"break" (whose probability is exponentially small any- 
way) should not modify the hopping energy window con- 
siderably. If the break is sufficiently far from the bound- 
ary it should not affect the current distribution among 
the boundary LSs. Hence the break mechanism has neg- 
ligible effect on hopping thermopower although it greatly 
modifies the hopping conductance)^ As a result, the vari- 
ance of the thermopower Var(e5r) = {{eST - {eST))'^) 
does not depend on the length of the system in the 
limit 2L — >■ oo. That is, it becomes a constant al- 
though (S) = for systems with particle-hole symme- 
try when 2L — > oo. In contrast, for the conductance, 
(InG) ^ -oo and Var(lnG) = ((InG- (InG))^) ^ 
when 2L cx))^ as the conductance is determined by 
the bulk and is significantly affected by the break mech- 
anism. 

We also computed the probability distribution func- 
tion of the thermopower for a specific set of parameters 
and plotted it in Fig. [7Jb). It is seen that the thcr- 
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mopower is mostly distributed in the range oi ^Eq < 
eST < Eq. The probability distribution function is not 
a normal distribution.— Rather, it has exponential tails, 
~ exp[—C\eS\T/Eo], at large \S\ with C being a constant 
depending on the parameters of the system. The expo- 
nential tails should come from the fact that LSs with high 
energies have exponentially small probabilities to be part 
of the backbone at the boundaries because the resistance 
between such sites and the lead is exponentially large. In 
the inset of Fig. Eljb) we also show how the thermopower 
evolves as a function of the system length for two random 
configurations. By increasing the system length there is 
an increasing probability to have, for example, a poorly 
conducting piece in the bulk of the system. It is seen that 
the thermopower saturates with large system length since 
the boundary effect is dominant. Meanwhile the different 
thermopowers for the two configurations vividly indicate 
the fluctuations of the thermopower even in very long 
systems. Relatively large mesoscopic (sample to sample) 
fluctuations in the low-temperature thermopower have 
been found also in the weak-disorder regimei^ 

B. Special contributions to the heat conductances 

In deriving Eqs. (|2T|) and ([22]) we have assumed that 
the heat current is carried by the percolating (span- 
ning) paths which transport the charge current as well. 
Therefore, the heat conductances K^, L3, and Kpe were 
all proportional to the conductance G. However, non- 
spanning paths can also contribute to heat conduction. 
This mechanism becomes especially important when the 
conductance G is sufficiently small. As the non-spanning 
paths do not conduct charge between the two leads, they 
have no contribution to G, Li, and L2. A non-spanning 
hopping path is schematically shown in Fig. [HKa). It is 
seen that by hopping back and forth between a lead and 
the nearby LSs having different energies, the associated 
phonon energy is transferred between the lead and the 
phonon bath. Therefore, whenever Tl ^ Tp {Tr ^ Tp) 
there will be heat fiowing between the left (right) lead 
and the phonon terminal. Even when there is no span- 
ning hopping path this scenario can lead to a finite heat 
conduction. Denoting the heat conductances due to non- 
spanning paths on the left and on the right sides hy Kl 
and Kr, respectively, the corresponding heat currents 
are Ql = Kl{Tp - Tl)/T and Qr = Kr{Tp - Tr)/T. 
Hence, when the non-spanning paths determine the heat 
conduction, we find using the definitions of currents and 
affinities in and below Eqs. [S] and the transport coeffi- 
cients of Eq. [SJ 

= (Kl + Kr)/4 , L3 = [Kb. - Kl)/2 , 
Kpe ^Kl + Kr = iK'^ . (33) 

In Fig. El^b) we show the numerically-computed aver- 
ages of K°/{GE^e-^) (curves with •), \L3\/{GE'^e-'^) 
(curves with A), and Kpe/{GE'^e~'^) (curves with □). 
It is seen that the heat conductances are too large to 



be explained by the contributions from the spanning 
paths, for which the upper bound is GE^e~^. This 
becomes more significant at lower temperatures or for 
longer lengths L, where the conductance is G is reduced 
but the non-spanning paths are only marginally affected. 
Clearly, the Wiedemann-Franz law totally breaks down 
here. These results confirm that the non-spanning paths 
can have much larger contributions to the heat conduc- 
tances than the spanning ones, when the conductance G 
is suppressed. Most of the non-spanning paths are also 
close to the boundary since the sites there are much more 
strongly coupled to the leads and the conductance and 
the number of long paths are substantially reduced)^ 
The boundary effect is demonstrated in Fig. [SJc) . One 
notes that the distance dependencies in this case are dif- 
ferent from those of the changes in S and Sp due to the 
different conduction mechanisms. Finally, this scenario 
is not important in ID NNH systems when the tunneling 
conductance between the left (right) lead and the LSs 
other than the leftmost (rightmost) LS is small enough. 
It can, however, play a role in 2D NNH systems when 
more than one LSs is strongly coupled to each of the 
leads. 



VI. CONCLUSIONS AND DISCUSSION 

The study of the thermopower in the hopping 
regime^iiS is augmented in this paper in two ways: 1. 
by considering the appropriate threc-termina l^^'^^ case; 
and 2. by studying the possibly-all-important effect of 
the edges on the thermoelectric transport. We empha- 
size that the three-terminal picture is dictated by the fact 
that hopping conductance necessitates energy-exchange 
with a (usually bosonic) thermal bath. 

We studied the boundary effect on the thermoelectric 
properties of finite ID and 2D hopping systems. We find 
that the boundary effect may play a crucial role for the 
thermopowers. This is first shown for a simple three-site 
hopping model and then for ID NNH systems via ana- 
lytical and numerical discussions. For ID VRH systems 
qualitative arguments and numerical results indicate that 
only the LSs with a distance from the boundaries smaller 
than or comparable to the Mott distance Rm can affect 
the thermoelectric properties considerably. As a conse- 
quence, the thermopowers of a specific sample of a very 
long ID hopping system where the particle-hole asym- 
metry is negligible on average is still finite and fluctuat- 
ing due to the insufficient averaging at the boundaries. 
This is confirmed by numerical calculations in ID VRH 
systems where a nonzero variance of the thermopower 
persists and eventually becomes a constant at very large 
system size. At the same time the average of the ther- 
mopower over many samples does vanish. We emphasize 
that the sample-dependent changes in the thermopower 
due to modifying the middle of the sample do exist, but 
they can be much smaller than the corresponding changes 
in the conductance. 
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For 2D systems we first found a situation which resem- 
bles the NNH ID cases: when the electronic leads are geo- 
metrically sharp and each of them is coupled with a single 
LS (as with a high- resolution STM probe). In this type 
of systems the boundary effect completely determines the 
thermoelectric properties. However, in other types of 2D 
hopping systems the bulk effect can also be important, 
but much less so than for the conductance. This is man- 
ifested by a simplified type of 2D hopping systems which 
are made up of parallel ID hopping chains where there 
is no hopping between different chains. The total ther- 
mopower is an average of the thermopower in each chain 
weighted by the conductance in that chain. Changing 
the central part of the system will alter the electric cur- 
rent in each chain differently. Although the thermopower 
in each chain docs not change, the weights of the vari- 
ous chains does change. Therefore the total thermopower 
is modified by changing the central part. This modifica- 
tion can again be much smaller than that for the conduc- 
tance (for example, when the latter is changed by many 
orders of magnitude by changing the bulk conductance 
that much) and it disappears upon ensemble averaging. 
In general the sample-specific thermoelectric properties 
depend mostly on the boundaries and in a limited fash- 
ion on the bulk whenever the current distribution at the 
boundaries can be affected by the bulk. This includes 
ID VRH systems, realistic 2D ones (for which some nu- 
merical results are presented) and we propose also 3D 
ones. 

For the thermal conductances sector of the linear- 
transport matrix [see Eq. ([6|)] for the three-terminal ge- 
ometry considered here we find new contributions. These 
are not due to the usual percolating paths, and will be- 
come important whenever the electrical conductance, due 
to the latter, is small enough. This implies a very serious 
breakdown of the Wiedemann-Franz law. 
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linear-response regime, 

_ \e,-e.,\ 5Tp 

Ul,bX£i) = \ TfT- ■ Al 

e el 

Here J/i; ^ fj,i ~ ^ and 5Ti ~ Ti — T with i ^ L, R, P. To 
simplify the calculation one may choose fj. — {hl +Mh)/2 
and T = {Tl + Tr)/2, so that SfiL = Si-Ir = Sfi/2, 
STl -STr = ST/2, and 6Tp = AT. According to the 
sign convention in Eq. ([9]), ±Uij = [sj — Si)dTp/eT. The 
final form of the equations to be solved is thus 



with 



An — ^ Gik + GiL + GiR, 
Aij = -Gy- (for i ^ j), 

' 5llL_ ^ - /i 5Tl 

e e T 

SfJ.R _^ Sj- ^J, STr 

e e T 

£i - £fc STp 



Zi — GiL 

+ GiR 



Eg- 



e T 



£j - M 



where 



li^L = GiL 
h^R — GiR^ 



SfJ.L 



flSTL 



e e T 

5hr ^ Ei - /i &Tr 
e e T 



(A2) 



(A3) 



Once the UiS are obtained by solving Eqs. (|A2I) . the 
three currents are found straightforwardly, 

i 



(A4) 



(A5) 



Appendix A: Numerical scheme for solving the 
resistor network in three-terminal geometries 

Here we present the numerical scheme for solving the 
resistor network in three-terminal geometries, in the 
linear-response regime. The key quantities to be calcu- 
lated are those appearing in Eq. Ui and Uij. In the 



The transport coefficients are obtained by computing the 
currents for three different cases: (i) Sfi but ST = 
AT = 0, (ii) ST ^0 but Sfi = AT = 0, and (in) AT ^ 
but (5/.J = ST = 0. Using Eqs. © in case (i) one 
obtains G, Li, and L2. In case (ii) one finds Li, K^, and 
L3, and case (iii) yields L2, L3, and -R'pe- The Onsager 
reciprocity relationships can then be verified from the 
numerical computation explicitly. 
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Appendix B: Comparing the conductance in the 
dominant hopping path with other conductances in 
the three-site NNH model 



Besides the dominant hopping path demonstrated in 
Fig. mja), there are the following transport processes: 
(A) elastic tunneling through the whole system; (B) 
tunneling from the left lead to LSI, hopping from LSI 
to LS2 and then tunneling into the right lead; (C) 
tunneling from the left lead to LS2, hopping to LS3 
and then tunneling into the right lead; (D) tunnel- 
ing from the left lead to LSI, hopping from LSI to 
LS3, and then tunneling into the right lead. One must 
keep in mind the assumption that Gil and Gsr are 
much larger than all other conductances. The conduc- 
tance of (B) and (C) are Gb = Gi2G2r/{Gi2 + Gzi?) 
and Gc = G2lG23./{G2l + G23). If G12 - G23 > 
G2ij; ~ G2L then such contributions can be negligible. 
According to Eq. pT|) this condition is fulfilled since 
G2L/G12 ^ exp[-2(r2L - ri2)/^] = exp(-2riL/0 « 1 
and G2fl/G23 ^ cxp[-2(r2fl-r23)/^] = exp(-2r3fl/0 < 
1. The conductance in process (D) can be considered 
similarly. 

It remains to consider the conductance of process ( A). 
The following analysis generalizes the one of Ref. [la 
which discusses only the two-LS assisted tunneling [see 
Eq. (jB3[) below]. The tunneling conduction consists of 
several contributions. It can be assisted by one, two, or 
three of the three LSs. For example, the tunneling con- 
ductance through LSi = (1, 2, 3) can be written as 



9i 



^\aiL?\aiR\^ PlPr 



(Bl) 



with PL and pR being the density of states in the left 
and right leads, respectively, jaiil ~ |ae| exp(— r^i/^) 
and \aiR\ ~ |ae| exp(— rjij/^) where \ae\ measures the 
tunnel coupling strength between the electronic states. 
The asymptotic behavior of gi is thus 



9i 



e^e^ '^\ae\'^PLPRCx.V 



4i 



(B2) 



where 2L = riL + riR is the length of the system. Simi- 
larly one can find the tunneling conductance through two 
LSs i^i = (1,2,3) as 



2 —2 — 2i ifi / 
m ^ e £i £j \ae\ PLPRCXp [ — 



(B3) 



Similar exponential dependence is also found for the tun- 
neling through three LSs. The asymptotic behavior of the 
total tunneling conductance Gtun is then 



G. 



e et„„|ae| PLPRexp 



AL 
T 



with e. 



(B4) 



In compar- 



ison, the conductance of the hopping channel is given by 



Eq. (Ull). If G12 - G23 and |ei - ^1, 1^2 - > hgT then 

-lep 



X exp 



|£i - mI + k2 - ^J'\ + |£i - £2! _ 2ri2 



2kBT 



It is seen that the hopping conductance is limited by the 
exponential factor at very low temperatures. Therefore 
the hopping conduction dominates at relatively high tem- 
peratures while the tunneling is more important at low 
ones. Ignoring the difference in the tunnel coupling and 
electron-phonon coupling, i.e. to a logarithmic accuracy, 
the crossover temperature between the two types of con- 
ductance, Tx, is given by 



ksTx 



H^i -^2 



4(2L - 7-12) 



(B5) 



Appendix C: A probability analysis of 
thermoelectric transport in the NNH three-site 
model 

Denoting the probability for an electron at LS 1 to be 
transferred to LS 2 per unit time by Pi^2, and that of the 
transfer from LS 2 to LS 3 by P2-S-3, the entire probability 
per unit time for an electron to be transferred from LS 1 
to LS 3 by passing LS 2 is 



Pi- 



(CI) 



P2^3 is the probability for the transfer from LS 2 to LS 
3 when the electron is already at LS 2. In that case 
there are two possibilities: the electron can either hop 
to LS 1 or to LS 3. Hence P2_j.3 is given by the ratio 

-F2^3/(-F2^i + ^2^3), and consequently 



^1^3 



Pi 



^2-^^2-5-3 



(C2) 



^2^-1 + ^2^3 

The probability per unit time for the reversed process is 



P3 



P3^2P2^1 



P2 



P2 



(C3) 



The Fermi golden-rule [see Eq. ([TJ] implies that Pi^2 = 

712/1 (1 - /2)iV21 , P2^1 = 712/2(1 - fl){N21 + 1), P2^3 = 

723/2(1 -/3)(A^23 + 1), and P3^2 =723/3(1 -/2)Af23. At 
equilibrium Pi->3 = ^3-5.1. When the system is out of 
equilibrium, in the linear-response regime, one has 

I = Pl-^3 — Pa^l 

Pl^2p2^a — Pz^2P2^1 



P2^1 
G12G23 r 



G12 + G23 

G12G23 r Sp 



P2^3 

Ui 



G12 ^~ G; 



23 



+ 



U12 + U3 + U23] 

^31 ST ^ UJ31 AT 



knT knT T k^T T 



(C4) 



This confirms the results obtained from the rate equation 
method, i.e., Eq. 
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FIG. 5. (Color online) ID VRH. (a): The boundary part of 
the hopping backbone, generalizing Fig. [2^. The gray region 
depicts the main part of the backbone. The connections 
(arrows) to the leads are through the LSs (red dots) near the 
boundaries, (b) and (c): Length {L — Rm)/Rm dependence 
of the average (| ln(|X'/X|)|). In (b): the conductance 
X = G (dashed curves) and the thermopower X = S (solid 
curves); In (c): the three-terminal thermopower X — Sp. 
Results given for three different temperatures ksT =15 (•), 
30 (A), and 60 (□). rg = 100, p = 0, ^ = 0.1, and p = 0.03, 
fcsTo = 333, Ec = 424, 600, and 849 for the three temper- 
atures respectively. Along each curve from left to right the 
number of LSs for the first, second, fifth data point are 
N = 30, 50, 70, 100, and 120 respectively. The Mott distances 
are Rm = 0.4,0.49,0.54,0.58,0.60(0.28,0.35,0.38,0.41,0.42) 
and the corresponding Mott energies are Em = 
120,147,161,174,180(170,208,227,246,254) for the five 
points at ksT = 15(30) respectively. Rm (Em) at fcsT — 60 
is half (two times) of that at ksT = 15 for the five points 
respectively. The average nearest neighbor distances 2L/N 
for the three T's are 0.039, 0.028, and 0.02 respectively. The 
results are averaged over 10® random configurations. The 
curves in (b) and (c) are guide to the eye. 
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FIG. 6. (Color online) The effect of decreasing the conduc- 
tances in the middle part of a random ID VRH system, by 
a factor of rg, on the conductance and the thermopower. 
X = G 01 S, with X'/X denoting the ratio of the con- 
ductance or thermopower after changing the middle, to the 
original value. The curves with A depict the thermopower 
while the ones with • depict the conductance. The parame- 
ters are: fcsT =15, = 0, ^ = 0.1, p = 0.03, ksTo = 333 
and Ec — 424. The average distance between adjacent LSs is 
0.039. For the longer (a) system (iV = 400, 2L = 22Rm) the 
resistances in the region { — 2Rm, 2Rm) are increased by a fac- 
tor of Tg. For the shorter (b) system (A'' = 40, 2L = 3.5i?M) 
the resistances in the region {—Km, Km) are increased in the 
same way. 




FIG. 7. (Color online) Thermopower inlD VRH. (a): The 
variance of the thermopower Var(e5T) = ((eST- {eST)f) 
as a function of the length of the system 2L/Rm at fcsT =15 
(•), 30 (A), and 60 (□). Along each curve from left to right 
the number of LSs for the first, second, seventh data 
point are N = 30, 50, 70, 100, 120, 200, and 400 respec- 
tively. The other parameters are the same as in Fig. [S] (b): 
Probability distribution function (PDF) of the thermopower 
for ID VRH systems, where ksT = 15, ^ = 0, C = 0.1, 
p = 0.03, Ec = 424, and A'^ = 120. The average distance is 
2L/N — 0.039. The results are obtained from 10® random 
configurations. The straight red line is an exponential fit to 
the tail, ~ exp( — C|eS'|T/i5o) with C ~ 5.5. In the inset we 
show how the thermopower evolves as a function of the sys- 
tem length for two random configurations. The parameters 
are: ksT =30, = 0, ^ = 0.1, p = 0.03, and Ec = 600. The 
average distance between adjacent LSs is 0.028. In figure (a) 
the curves are guide to the eye. 
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FIG. 8. (Color online) (a): Schematic of the non-spanning 
hopping paths which may contribute significantly to the heat 
conductances. The wavy lines denote the phonons involved 
in the processes, (b) and (c): Distance {L — Rm)/Rm depen- 
dence of the average of (b) \X\/{GE^e-^) with X = K° (•), 
X = L3 (A), and X = Kp^ (□) for ksT =15 (solid curves) 
and 30 (dashed curves), (c) |ln(|X7X|)| for X = Kp^ (solid 
curves), X = (dashed curves), and X = L3 (dot-dashed 
curves) for three different temperatures with fesT =15 (•), 
30 (A), and 60 (□). The other parameters are the same as 
those in Fig. [S] The results are obtained by averaging over 
10® random configurations. Note that in figure (c) the curves 
for K° are quite indistinguishable from those of Kp^ since 
Kpe = 4_R'g when the non-spanning paths determine the heat 
conduction. The curves in figures (b) and (c) are guide to the 
eye. 



